Effect of deconfinement on resonant transport in quantum wires 
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The effect of deconfinement due to finite band offsets on transport through quantum wires with 
two constrictions is investigated. It is shown that the increase in resonance linewidth becomes 
increasingly important as the size is reduced and ultimately places an upper limit on the energy 
(temperature) scale for which resonances may be observed. 
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I. INTRODUCTION 

Recent technological advances have enabled semicon- 
ductor nanostructures to be fabricated with feature sizes 
down to tens of angstroms. Such structures include ar- 
rays of 'self-organized' quantum dots and quantum wires. 
The former are grown by heteroepitaxial deposition in 
which the over-layer material has the larger lattice pa- 
rameter and forms quantum dots to relieve clastic strain 
00. The latter may be achieved by heteroepitaxial 
growth on 'v'-groove surfaces, produced by optical lithog- 
raphy and etching, for which the over-layer atoms dif- 
fuse preferentially towards the base of groove producing a 
quantum wire with a crescent-shaped cross section [|3| [10| . 
These structures have potential opto-electronic applica- 
tions, such as light-emitting diodes, low-threshold lasers 
and single-electron devices. 

In this paper we consider the ballistic transport of elec- 
trons through quantum wires in which there are two con- 
strictions defining a small region of quasi-confincmcnt. 
Such a system can have a rich and complex resonance 
structure |ll]-[l5|]. It has been investigated in detail by 
Nakazato and Blaikie |pd| , who consider characteristic 
sizes down to ^100 nm and assume that the quantum 
wire is defined by infinitely high barriers. Such an as- 
sumption is reasonable provided that the resonance lev- 
els considered are at sufficiently low-energy, i.e. some- 
what lower than that of the real band edge in wider-gap 
material of the heterojunction. For the wire geometry 
they consider, the resonances are very sharp at T = 0. 
Inter-channel (mode) mixing of the electron waves is im- 
portant for higher-order resonances and can give rise to 
unusual effects such as anti-resonances. At these length- 
scales, however, the possible observation of such reso- 
nances would be restricted to very low temperatures, 
where they would be easily swamped by defects and dis- 
order. 

The main purpose of the present work is to investi- 
gate the resonance structure for very narrow wires for 



which the energy (temperature) scale is higher and may 
even approach room temperature. Such a possibility is, 
of course, extremely important for device applications. 
As the size shrinks, the effect of a finite barrier for an 
electron in the wire becomes increasingly relevant, the 
main effect being to broaden the resonances and restrict 
their number. Indeed, for a sufficiently small confining re- 
gion, only a single isolated resonance level remains, with 
a continuum at higher energies. The increasing width 
of the resonance level is due to deconfinement of the 
electron wave-function into the classically forbidden re- 
gion beyond the potential step. However, all energies 
increase with increasing confinement and the important 
criterion is whether the lowest resonance peak can be re- 
solved from the next resonance peak (or the continuum), 
which sets the temperature scale. To be specific we shall, 
in what follows, mainly consider two-dimensional struc- 
tures with a conduction-band offset of 0.4 eV, which is ap- 
proximately the maximum offset in the conduction band 
(direct gap) for GaAs wires on an Al K Gai_a; As. 



II. MODEL AND METHOD 

For simplicity, we restricted the model to two dimen- 
sions with confinement in the y-direction and propa- 
gating in the a;-direction. The wire shapes under con- 
sideration are symmetric around the x axis and shown 
in Fig. 1. The width of the wire is parametrized as 
a(x) — ao — ai sin 2 2ttx/ ' ai for < x < a,2 and a(x) = ciq 
otherwise. For comparison we choose two generic geome- 
tries. In Fig. 1(a) the confined region is only weakly 
coupled to the 'leads' and we thus expect strong (nar- 
row) resonances. This structure is expected to behave in 
a similar way to a quantum dot which is separated from 
its leads by tunnel barriers. On the other hand, for the 
geometry of Fig. 1(b), the constrictions are much smaller 
resembling a wire with weak thickness fluctuations. 

We model further the wire as two regions of constant 
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potentials, V = within the wire 'boundary', and con- 
fining potential V = Vq > outside the wire. The corre- 
sponding two-dimensional Schrodinger equation reads 



n 2 ( d 2 



u 2 



2m* \dx 2 dy 2 



*(a?, y) + V9(x, y) = EV{x, y). 



(1) 



Here the electron effective mass m* is chosen to be that 
of GaAs, equal to 0.067 times the free electron mass and 
E is electron energy relative to the conduction band edge 
in the wire. 




20 nm 



FIG. 1. The width of the wire is always parametrized as 
a(x) — ao — ai sin 2 l-Kxjai. (a) Weaker coupling: ai = 0.8ao 
and a,2 = 2ao- (b) Stronger coupling: ai = 0.4ao and 
a,2 = 2ao. (c) Shape of the 'dot' fixed with a± — 5 nm, 
a,2 = 20 nm and various wire thickness ao = 6 — 10 nm. 

The wave function is expanded in elementary modes 
or channels pr| , 



(2) 



where the basis functions <& n (x,y) are orthogonal solu- 
tions of one-dimensional Schrodinger equations in the y 
direction for fixed x, with eigenenergies E n (x). 

Substituting the expansion Eq. (|2|) into the Schrodin- 
ger equation Eq. ([!]) leads to a set of ordinary linear dif- 
ferential equations for the channel n wave function ^ n (x), 
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where k 2 n (x) = 2m* [E ~ E n {x)]/h 2 and 



d 2 $ m (x,y) 
dx 2 



dy, 



b nm (x) = -2 J ^ n (x lV ) d ^ v) dy. (4) 

For the geometry considered here, there are no real bound 
states in the x -direction; only scattering states are rele- 
vant. 

Zero temperature conductance is calculated using the 
Landauer formula [fPi 



G = Go T(E), 



(5) 



where Go = 2e 2 /h and E is here the Fermi energy of 
the electrons in the leads. The transmission probabil- 
ity, T(E) , is the sum of transmission probabilities for all 
channels, n, at energy E, i.e., 



T(E) = "£T n (E), 



where 



T n (E) 



\tnm{E)\ 



(6) 



(7) 



and t nm {E) is transition amplitude for scattering from 
channel n to channel m. 

At finite temperatures the conductance is calculated 
using a generalized Landauer formula [fl9[ 



G — Gr 



df(e-E,T) 
de 



de, 



(8) 



where f(e,T) = [1 + exp(e//cBT)] _1 is the usual Fermi 
function. This form is also based on the assumption that 
motion within the wire is ballistic, the effect of temper- 
ature being merely to change the energy distribution in 
the leads, thus allowing electrons above and below the 
Fermi energy to contribute to the conductance. For nar- 
row wires, the tail of the Fermi distribution with energy 
e > Vq can be significant and the contribution of these 
electrons to conductance will depend on the size and 
properties of the barrier region. If this region is large, 
then electron motion will be diffusive and may be de- 
scribed by an effective conductivity, proportional to 5n, 
the number of electrons in the Fermi tail. The conduc- 
tance due to these electrons will be further inhibited by 
rough-surface scattering for mesa structures produced by 
etching. In this paper we shall only consider the conduc- 
tance due to electrons within the wire by introducing an 
energy cut-off at e = Vq, as shown in Eq. (Sh. Thermally 
activated electrons in the barrier region, which give rise 
to a series conductance, will be considered in future work. 

As pointed out in Ref. O], the Schrodinger equation 
Eq. (0) and the expression for conductance, Eq. (|| ), are 
invariant under the scaling transformation, 
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x, y -> Xx, Ay, 
E,V ^ \- 2 E,\- 2 V, (9) 
T -> A~ 2 T. 

To solve the system of differential equations, Eq. (^), 
we first fix the number of channels, N, and then de- 
termine the eigenfunctions <& n (x, y) with corresponding 
eigenenergies E n (x) for n < N . N must, of course, be 
sufficiently large to ensure convergence. For narrow wire 
this poses a problem since channels with energy above 
the barrier, Vq, form a quasi-continuum and because of 
inter-channel coupling [c.f. terms at the r.h.s. of Eq. (||)], 
these high-energy channels can have a significant effect 
on the eigenstates of electrons confined to the wire. This 
may be understood in a perturbation theoretic sense. 
An electron in state ip n {x) with E n (x) < Vq may make 
a virtual transition to a state ip m {x) in the continuum 
(E m > Vq) and such transitions become very important 
for Vo — E n small, i.e. for confined states close to the top 
of the barrier. These excursions into the quasi-continuum 
enhance deconfinement into the classically forbidden re- 
gions. This is particularly important near narrow con- 
strictions, such as the 'necks' in figure 1(a). The main 
effect of this, and the leakage of the base states ip n into 
the barrier region, is to broaden resonances compared 
to cases with Vo = oo, considered in Ref. To 
model the quasi-continuum we introduce infinite barriers 
at y = ±L/2. The basis functions, <fr n (x, y), x-fixed, are 
the simple standing waves. Some care is needed in opti- 
mizing N and L to ensure convergence. L must be made 
sufficiently large to ensure that the effect of finite barriers 
on the confined wave-functions is negligible. On the other 
hand, L cannot be made too large otherwise the number 
of required channels with E > Vq becomes impractical. 
For the cases we considered, L < 20ao and N < 30 were 
sufficient to ensure convergence. The analytic expres- 
sions for the functions & n (x, y) enable the coefficients 
a n m(x) and b nm (x) in Eq. (J3j) to be computed efficiently. 
The system of differential equations Eq. (H|) belongs to 
a class of 'stiff' equations for which direct integration is 
generally not stable. The basic problem is exponentially 
increasing solutions with imaginary wave vectors of dif- 
ferent orders of magnitude leading to round off errors 
and divergent results. This problem was solved for the 
present system by dividing the wire along ir-direction in 
M sections. For each section we first determined 2N in- 
dependent solutions using the fifth order Runge-Kutta 
numerical method. The length of each section was cho- 
sen to ensure stable numerics in that section, the limit- 
ing factor being the number of channels with imaginary 
k n (x). In our case up to 10 such channels were taken 
into account within each section, with M < 10 sections. 
Matching the solutions at each boundary yields sets of 
linear equations from which the transition amplitudes, 
tnm(E), may be determined. 



III. RESULTS 

Fig. 2(a) shows electron conductance at T — versus 
energy for the dot geometry of Fig. 1(a), i.e., a\ = 0.8ao 
and ai = 2eto for wire widths ao from 4 nm to 20 nm and 
the barrier height Vo = 0.4 eV. The open circles, o, in this 
and the remaining figures correspond to a Fermi energy 
E = Vo, as explained in the previous section, electrons 
with energy greater than this are not 'bound' to the wire 
and their conductance will be dominated by the prop- 
erties and size of the barrier region. For convenience of 
presentation and to emphasize the effects of scaled units 
of energy, we choose scaled units of energy, E/E , where 
E = 7i 2 /(2m*a§), the ground-state energy for an elec- 
tron in a one-dimensional well of width ao with infinitely 
high potential walls. For perfect confinement (Vo = oo), 
the scaling invariance, Eq. (^|), at T = gives, 

G(Aa , E/E (\a )) = G(a , E/ E {a )) (10) 

and hence all wires of the same shape have identical con- 
ductance curves for Vo = oo 111]]. This universal curve is 
approximately that for ao = 20 nm in Fig. 2(a) (for which 
the finite barrier height is irrelevant). We note that the 
(first) resonance occurs at E <~ 2.5-E'o, which may be in- 
terpreted approximately as a resonant bound state with 
energy Eq due to confinement in the y-direction and en- 
ergy 1.5i?o due to quasi-confinement in the x-direction. 

We see that the effect of finite Vq is to shift the position 
of the first resonance to lower energies with decreasing 
wire width. This is, of course, due directly to the decon- 
finement effect of a finite band offset. For very narrow 
wires the resonant level is pushed towards the continuum 
(i.e. unbound in 2D) at energy E = Vq = 0.4 eV and only 
one resonant level is possible. This is seen to be the case 
for wires of width 4, 5 and 7 nm. The corresponding en- 
ergies from the peak to the continuum are 83, 151 and 
232 meV for these wires. For wider wires the energy split- 
ting between the first two peaks reduce in energy, eventu- 
ally decreasing like a^ 2 (approximately). Furthermore, 
the resonances always become broader with decreasing 
wire width, again a consequence of the deconfinement 
effect of finite band offset. Thus, at zero temperature, 
multiple resonances with the highest resolution (ratio of 
resonance separation to resonance width) occur for wide 
wires, the limiting resonance width being determined by 
the geometry (strength of the effective tunnel barriers 
into the confined region). However, the absolute energy 
scale is, of course, very small and these sharp resonances 
are rapidly broadened with temperature. This is shown 
in Fig. 2(b) where for the same wires we plot conductance 
at T = 100 K and 300 K with thick and thin lines, respec- 
tively. Open circles o again represent an energy cut-off 
E = Vo, above which conduction is primarily through the 
barrier region. The thermal broadening is accompanied 
by a reduction in the peak heights which are barely re- 
solvable at room temperature for wires with oq > 10 nm, 
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though clearly resolvable and optimum for ciq ~ 5 nm. 
This is significantly better than the 4 nm width wire for 
which the proximity of the continuum has a large effect. 



1.5 - 4nm(0.4eV) " 

" (a) T=0 -—5(0.63) 
7(1.23) 

; — 10(2.5) ; 




in the third dimension. The behavior of such wires at 
low-temperatures, i.e. conductance versus energy/gate 
voltage for various widths, would be the same as that of 
narrow wires at higher temperatures. |^] 

In Fig. 3 we show conductance plots at absolute zero 
and T = 100 K for the more weakly confining wire of 
Fig. 1(b), for which a\ = 0.4a and 02 = 2ao- The 
behavior is seen to be qualitatively similar to the more 
strongly confined geometry wire though the resonances 
are broader, reflecting the weaker confinement along the 
wire (smaller effective tunnelling barriers). However the 
thermal broadening is still largely governed by the overall 
wire thickness. For example, if we compare the 5 nm 
wires in Fig. 2(b) and Fig. 3 at T = 100 K (thin lines) 
we see a similar relative increase in the halfwidth and 
decrease in the peak height compared with absolute zero. 
This shows that even for small thickness fluctuations, 
resonance peaks can persist to quite high temperatures 
for narrow wires. 



1.2 
1.0 
0.8 
^0.6 
0.4 
0.2 
0.0 



1 1 

- (b) 


1 1 1 1 1 

4 nm 


'_ j\\ — - 10 




1 \ 

\ 1 \ - 
\ 1 1 




V / 1 

/ \\\ 1 ,' ! 

\ y y / - 





1.5 



1.0 



55 



0.5 



0.0 



12 3 4 



FIG. 2. (a) T — conductivity [geometry from Fig. 1(a)] 
for various ao — 4 — 20 nm and fixed Va = 0.4 eV (or, equiva- 
lently, with fixed ao = 4 nm and various Vb = 0.4—10 eV). (b) 
Conductivity at T = 100 K and 300 K presented with thick 
and thin lines, respectively, for fixed Vb = 0.4 eV and various 
ao = 4 — 10 nm. Circles o represent upper limit E = Vb of 
calculation. 

Because of the invariance under scaling, Eqs. (|J), the 
same behavior occurs for larger wires with a smaller band 
offset, though the overall energy scale is lower. Thus 
Fig. 2(a) may be regarded as a plot of conductance for 
wires with same width (4 nm in this case) but differ- 
ent band offsets, Vb- Increasing Vq increases the con- 
finement and has the same effect as increasing the wire 
width with fixed Vb and visa versa, apart from an overall 
change in energy/temperature scale. Indeed, the effects 
of deconfinement could be investigated experimentally 
for relatively large wires by fabricating quasi-2D wires 
with small conduction-band offsets. The quasi-2D be- 
havior would be achieved by ensuring high confinement 
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FIG. 3. T — conductivity [geometry from Fig. 1(b)] for 
various ao = 5 — 1000 nm and fixed Vb = 0.4 eV is pre- 
sented with thick lines. Thin lines represent the correspond- 
ing T = 100 K result. Circles o represent upper limit E — Vo 
of calculation. 

Finally we consider the effect of reducing only the 
thickness of a quantum wire, ao, whilst otherwise main- 
taining the same shape, a\ = 5 nm and 02 = 20 nm, as 
shown in Fig. 1(c). This has the effect of producing a 
quasi-lD quantum dot. There are two competing effects 
as the wires is made narrower: the effective tunnel 



* This assumes very clean wires and parabolic bands. In 
practice one would have to take into account non-parabolicity 
effects for very narrow wires and the effects of disorder would 
become increasingly important for larger wires. 
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FIG. 4. (a) T = conductivity [geometry from Fig. 1(c)] 
for various ao = 6 — 10 nm and fixed Vb = 0.4 eV. (b) Con- 
ductivity at T — 100 K and the same geometry. Circles o 
represent upper limit E = Vo of calculation. 

barriers increase and there is enhanced deconfinement 
near the necks. It turns out that the increase in effective 
tunnel barriers is the greater effect, resulting in sharper 
resonances as the wire thickness is reduced, as shown 
in Fig. 4(a). This should be contrasted with Fig. 2(a), 
which always gives rise to a resonance broadening when 
the overall size is reduced. However, we point out that, 
unlike the case of thick wires, this reduction in linewidth 
is significantly less than that given by a single-channel 
approximation. This is because there is coupling to the 
2D continuum near the neck region. This coupling de- 
pends on both the rate of narrowing of the wire and its 
absolute width. The other main effect shown in Fig. 4(a) 
is the appearance of further resonance peaks below the 
continuum. This reflects the quasi- ID nature of the con- 
finement region, the higher-lying peaks corresponding to 
higher harmonics along the length of the effective quan- 
tum dot. Their separation is set by the length of the 
confinement region along the wire. Conductivity at fi- 
nite a temperature T = 100 K is plotted in Fig. 4(b). We 



see that for the wire thicknesses considered, the higher 
harmonics broaden rapidly with increasing temperature, 
merging into the quasi-2D continuum. However, it can 
be seen that the lowest resonance remains distinct. 

IV. CONCLUSIONS 

With realistic conduction band offsets, the lowest reso- 
nance peak for ballistic transport through quantum wires 
of fluctuating thickness giving rise to a quasi quantum 
dot should be discernible around room temperature for 
smallest structures close to the limits of present fabrica- 
tion techniques. Nevertheless, deconfinement effects due 
to finite band offsets are significant and ultimately the 
limiting factor for sufficiently small structures. In prin- 
ciple the resonances could be made to survive to higher 
temperatures if heterojunctions could be fabricated with 
even larger band offsets. The restriction to 2D in the sim- 
ulations is somewhat artificial here though the extension 
to true 3D with circular cross-section wires is feasible and 
calculations are in progress to estimate the expected en- 
hanced deconfinement effect which they would produce. 
However, the general behavior is expected to be simi- 
lar to that described in this paper. Other geometries, 
such as might be produced from a self-organized quan- 
tum dot connected to source and drain contacts 'ver- 
tically' is also expected to behave in a similar fashion, 
though there would be quantitative differences of course. 
Future work will consider other effects which become in- 
creasingly important for very small structures, including 
non-parabolicity and Coulomb blockade, particularly the 
effect of deconfinement on the latter. 
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